The classification of real singularities using Singular 
Part I: Splitting Lemma and Simple Singularities 



Magdaleen S. Marais and Andreas Steenpafi 

Abstract. We present algorithms to classify isolated hypersurface singu- 
larities over the rea l numbers according to the classification by V.I. Arnold 
l lArnold et al.l.ll985^ . This first part covers the splitting lemma and the sim- 
ple singularities; a second part will be devoted to the unimodal singularities 
up to corank 2. All algorithms are implemen ted in the Singular library 
realclassify.lib l lMarais and Steenpal . l2012h . 



1. Introduction 



Arnold et all (|l985l ) present classification theorems for singularities over the 
complex numbers up to modality 2 and for singularities over the real numbers up 
to modality 1, including complete sets of normal forms. For the complex case, he 
also gives an algorithm how the type of a given singularity can be computed, called 
the "determinator of singularities" (cf. [Arnold et al. . 19851 ch. 16), but this question 
is left open for the real case. The goal of this paper, together with its second part, is 
to present algorithms for the classification of isolated hypersurface singularities up 
to modality 1 and corank 2 over the real numbers with respect to right equivalence. 

We consider real functions with a critical point at the origin and critical value 0, 
i.e. functions in m^, where m denotes the ideal of function germs vanishing at 
the origin. Two function germs /, 5 G C M[[a;i, . . . , a;„]] are considered as 
right equivalent, denoted by f'^g, if there exists an M-algebra automorphism cj) 
of R[[xi , . . . , Xn]] such that 0(/) — g. 

We have implemente d all the algorithm s presented here in the computer al- 
gebra system SINGULAR ( Decker et al. . 20121) . The implementation is freely avail- 
able as a Singular library called realclassify.lib which relies on Singular's 
classify. lib to determine, for a given polynomial, the type in Arnold's classifi- 
cation over the complex numbers. The methods used in cla ssify .lib wil l not be 
discussed in this paper. For more information in this regard. iKriigerl (jlQQTI ) can be 
studied. 

In section ^ we introduce basic invariants and methods which are frequently 
used for the algorithmic classification in the subsequent sections. We first recall 
some basic results on the Milnor number and the determinacy in subsections 12.11 
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and 12.21 and we also recall how these invariants can be computed. As a further 
prerequisite, we show that the homogeneous parts of lowest degree of two right 
equivalent functions factorize in the same way over R (section 12.31 Proposition [5]) . 
We also show that in some cases, this factorization can even be carried out over Q 
which is important for the algorithmic aspect (Lemma [7]) . 

Using the Splitting Lemma (Theorem [9)) , any function germ / over the real 
numbers with an isolated singularity at the origin can be written, after choosing a 
suitable coordinate system, as the sum of two functions of which the variables are 
disjoint. One of the functions, called the nondegenerate part of /, is a nondegenerate 
quadratic form and the other function, called the residual part of /, is an element 
of m'^. The number of variables in the residual part is equal to the corank of 
/, denoted by corank(/). In this paper we only consider germs with corank 0, 1 
and 2. A version of the Splitting Lemma for singularities over M and a corresponding 
algorit hm are discussed in section |31 



In Arnold et al 



the real singularities of modality and 1 are classified 
up to stable equivalence into main types which split up into more subtypes depend- 
ing on the sign of certain terms. Two functions are stably equivalent if they are right 
equivalent after the direct addition of nondegenerate quadratic forms. Hence after 
applying the Splitting Lemma, we only need to consider the residual part in order 
to compute the correct subtype. It can be easily seen that the subtypes are complex 
equivalent to a complex singularity type of the same name as its corresponding real 
main singularity type (see Table[T]). In fact there is a bijection between the complex 
types of modality and 1 and the main real types. Thus, if we can determine the 
complex type of a function germ, we only need to determine the correct subtype of 
the corresponding main real type. The classification of the residual part is given in 
section 131 together with explicit algorithms for each singularity type. 

2. Prerequisites 

2.1. The Milnor Number. We briefly recall the following well-known defi- 
nition: 



and p S Ag, the Milnor number of f at 



fi{f, p) := dimM \ M\[xi - pi, . . . , x„ 



'Pn\ 



e NU{oo} . 



Definition 1. For f e 

p is defined as 

df_ df_ 

dxi ' ' dxn , 

If p is the origin, we simply write fi{f) instead of ii{f,p). 

The Milnor number is known to be finite at isolated singularities (cf. lGreuel et al.l 

20071 Chapter I. Lemma 2.3) and to be invariant under right equivalence (cf. 
Lemma 2.10 ibid.). It is thus an important tool for the classification of isolated 



singularities. We refer to lGreuel et al.l (|2007[ ) for more properties of this invariant. 

There is a well-known algorithm f or the computation of the Milnor number 
which is implemented in SINGULAR, see iGreuel and Pfisteij (|2008f ). pp. 526-528. 



2.2. The Determinacy. In general, the singularities we deal with in this 
paper are defined by power series, but algorithmically, we want to work with poly- 
nomials. It is thus important for our algorithmic approach that any power series 
defining an isolated singularity is right equivalent to a polynomial which can be 
obtained from it by leaving out terms of sufficiently high order. 
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Definition 2. Let f G M[[a:i, . . . , be a power series. 

(1) Let f = X^j^o fj decomposition of f into homogeneous parts fj of 
degree j . For k £N, we define the fc-jet of f as 

k 

jet(/,fc) :=^/.. 

1=0 

In other words, the k-jet of f can be obtained from f by leaving out all 
terms of order higher than k. 

(2) / is called fc-determined if 

Vgem^'+i; / ^ jct(/, fc) + 5 . 

The determinacy is, just as the Mihior number, both invariant under right 
equivalence and fin ite for isolated singularities. We cite the following statement (cf. 
Greuel et all I2OO7I Chapter I, Supplement to Theorem 2 . 23) du e to its importance 



for the algorithmic approach and refer to iGreuel et al. ( 2007 ) for further results 
regarding the determinacy: 



PROPOSITION 3. Let f e m cR[[xi, ... , a;„]] . // 

'9/ dj 

dxi ' ' dx 



Xi,...,Xn\] 



holds, then f is k-determined. 

As a consequence of this, any power series / which has an iso lated singu- 
larity at the origin is {n{f) + l)-determined (cf. iGreuel et all l2007l Chapter I, 



Corollary 2.24). But we can often compute a much better upper bound for the 
determinacy by using the above statement as in Algorithm [TJ 

Algorithm 1 Determinacy 

Input: / G . . . , x„] with an isolated singularity at the origin 

Output: an upper bound for the determinacy of / 
1: k := Milnor(/) + 1 

3: compute a standard basis G of (m^ J) w.r.t. a local monomial ordering < 
4: for (/ = 1, . . . , fc - 1) do 
5: if (NF< (m'+i , G) = 0) then 
k / 
break 
return k 



Remark 4. In Algorithm[l] the for-loop computes the minimal k gN such that 
the condition in Proposition [3] holds. Th is number is equal t o the degree of the so- 
called highest corner of (G) = (m^ J) (cf. lCreuel and Pfisteil . [20081 Corollary A.9.7) 
and can thus also be computed by combinatorial means with the Singular com- 
mand highcornerO which is often much faster. 

It is worth to note that both the Milnor number and the determinacy of an 
arbitrary power series / G K[[a;i , . . . , Xn]] do not change if we regard / as an element 
of C[[a;i, . . . jXra]]. The same holds for the output of the corresponding algorithms 
presented here. 
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2.3. Results regarding the factorization of homogeneous polynomials 
over M and Q. 

Definition 5. Let <j) he an M-algebra automorphism o/ . . . ,x„]]. For 
j > we define the j-jet of (j), denoted by (j)j , to be the automorphism given by 

<j>j{xi) -.^ iet{<j>{x^),j) \/i^l,...,n. 

The next result is in many cases a starting point for the algorithmic classifica- 
tion of the residual part, see sectional Given / and g with f g, it can be used to 
determine (f>i for some automorphism (j) such that (!){/) — g. 

Proposition 6. Let f,g£ M.[[xi, . . . be two power series with f g and 
k := ord(/) > 1. Let (p be an W-algebra automorphism ofW[[xi, . . . such that 

//jet(/, fc) factorizes as 

jet(/,A:) = /f •••/r 
inM.[xi, . . . ,Xn], then iet {g,k) factorizes as 

iet{g,k)^Mhr---Mftr- 

Proof. By assumption we have that f = fl^ ■ ■ ■ ft'' + /', where fl^ ■ ■ ■ f^* is 
homogeneous of degree k and the order of /' is greater than k. We denote the 
higher order parts of by 0* :— 4)~ 4>i. Since </) is a homomorphism, it follows that 

</>(/) = '^(/r---/r) + 0(/') 

= •••/*') +'/'(/') 

where (piifi^ ■ ■ ■ ft') is homogeneous of degree k and both (p*{fi^ ■ ■ ■ ft') and (p{f') 
are of order higher than fc. Hence 

jet(<?, fc) = jetW), fc) = Mfi' ■■■ft')- MfiY' ■ ■ ■ Mft)" ■ 

□ 

Since we do not want to work with rounding errors nor field extensions in the 
implementation of the proposed algorithms, the above result would not be of much 
help for this purpose without the following result. 

Lemma 7. /// G Q[x,y] is homogeneous and factorizes as 

(i) gf or (ii) gig^, 

where 51,52 G K[a;, y] are polynomials of degree 1 and d > 1, then f factorizes as 

(i) ag[^ or (ii) ag[g'^ , 
respectively, where g'i,g2 ^ Qt^^i?/] i^^e polynomials of degree 1 and a £ Q. 

Proof, (i) Let / = [aix + a2y)'^, 01,02 G M. Without loss of generality, 
suppose ai ^ 0. Then f = af{x + ^y)'^- Since the coefficient of a;** in / G Q[x, y] is 
af, we have af £ Q and therefore {x+ ^y)'^ G Q[a;, j/] which, by dehomogenization, 
leads to [x + ^)'^ e <Q[x\. Since Q is a perfect field it follows that G Q. Thus 
/ = ag'i, where a := af G Q and g'^ ^ x + ^y G Q[x, y]. 

(ii) Let / {aix + a2y){a3X + a^y)'^, ai, . . . ,04 G M. Suppose 01,03 ^ 0. For 
the cases 01,04 ^ 0, 02,03 7^ and 02,04 ^ the proofs are similar. We have 
O1O3 G Q analogously to part (i). Hence (x + ^y){x + f|j/)'* G Q{x,y] which in 
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turn implies {x + + G Since Q is a perfect field it follows that the 

roots of this polynomial are rational. Therefore / = cLg'ig'-z with a := aiaf e Q, 
g[ := (x+^y)eQ[x,yl and := {x + l^y) e Q[x,y]. □ 

3. The Splitting Lemma 

Definition 8. For f G M[[a;i, . . . ,a;„]], we define the corank of f, denoted by 
corank(/), as the corank of the Hessian matrix H{f) at 0, i.e. 

corank(/) :— corank(i/(/)(0)) . 

The following well-known theorem, called the Splitting Lemma, allows us to 
reduce the classification to germs of full corank or, algorithmically, to a polynomial 
contained in m'^ nR[a;i, . . . , Xc] for a given input polynomial of corank c. We present 
a version for singularities over the real numbers, taking into account the signs of 
the squares. 

Theorem 9. // / G rri'^ C ^[[xi, . . . ,x„]] has an isolated singularity and if its 
corank is c, then 

c+A n 

f^g- y: -1+ E 

i=c+l i=c+A+l 

with g G nM[[xi, . . . ,a;c]]. g is called the residual part of f and X is called 
the inertia index of f. Both A and the right-equivalence class of g are uniquely 
determined by f. 

The foll owing proof is based upon the proofs of the Theorems 2.46 and 2.47 in 
Chapter I of iGreuel etHI (|2007t) . 

Proof. The corank of the Hessian matrix of / at is c, so by the theory of 
quadratic forms over R there is a transformation matrix T such that 

• i?(/)(0) • T = diag(0, . . . , 0, -1, . . . , -1, 1, . . . , 1) . 
Therefore the linear coordinate change (a:i, . . . , x„) i— > (xi, . . . , x„) • T* transforms 
the 2-jet of / into (^— J2i=c+i ^1 + T^l=c+\+i where A is the inertia index of /. 
Applied to /, this transformation leads to 

/(3)(xi,...,x„) :=/((a;i,...,x„)-T*) 

c+A 

= 53 - 

2— c+1 1— c+A+1 i—c+1 

with ^3 G nM[[a:i, . . . , Xc]] and h^^^ G m^. The coordinate change (^^^^ defined 

by 

\xi-\hf\ i>c. 

yields 

/W(xi,...,x„):=/(3)(^(3)(xi,...,x„)) 

C+A n n 

i— c+1 2— c+A+1 i— c+1 



E E E -^-^^^ 
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with g4 G m"* nM[[a;i, . . . ,Xc]] and G m'^. Continuing in the same manner, the 
last sum will be of arbitrarily high order. It can be eventually left out because / is 
finitely determined as an isolated singularity. □ 

Since this proof is constructive, we can immediately derive Algorithm [2] from 

it. 

Algorithm 2 Algorithm for the Splitting Lemma 

,2 



Input: feme Q[xi, . . . , x„] and k such that / is fc-determined 



Output: the corank c of /. the inertia index A of / and g G n Q[a;i, . . . , ccc 



such that 

c+A n 

„2 



■i— c+l 2— c+A+1 

1: compute a transformation matrix T G Q"^^^ such that 

T* • H{f){0) ■ T = diag(0, . . . , 0, -1, . . . , -1, 1, . . . , 1) iV 

2: c := number of zeroes on the diagonal of 

A := number of entries equal to —1 on the diagonal of N 
P\x,,...,xn) :-/(r-(a;i,...,x„)*) 
for {I = 3, . . . , k) do 
write /^'^ as 

/ c+A n n 

/^'^-E5.- E E E -^-f^' 

j—3 i^c+l i— c+A+1 i—c+1 

with gj e n Q[xi, . . . ,Xc] and hf^ e m'^-^ 
7: (yS^'H/^'^) where ^('^ is defined by 

^(;)^^ -J I ^^i, i < c, 

\xi - lhf\ i > c. 

3=3 9} 

9: return c, A, g 



4. The real classification of the residual part 

In Arnold et al. I (|l985h . the real singularities of modality are divided, using 
stable equivalence, into main types which split up into more subtypes by changing 
the sign in front of certain terms. It is known t hat the modality does not decrease 
under complexification (cf. [Arnold et all . Il985l pp. 273-274). So by applying the 
algorithms for the complex classification to the real normal forms, it is easy to see 
that in modality 0, there is a one-to-one correspondence between the complex types 
and the real main types. As we will see in the second part of this paper, the same 
holds true also in modality 1, but in both cases, this is not clear a priori and can 
only be deduced from the independently derived complex and real classifications. 
In fact, it is not known whether the modality is preserved under complexification 
in general. 
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Both the real and complex normal forms of the simple singularities are listed 
in Table [TJ From here onwards we will work with stable equivalence. For all 
degenerate forms it is thus only necessary, after applying the Splitting Lemma, to 
consider their residual parts, i.e. germs in m'^. 

Table 1. Real normal forms of singularities of modality 0. 





Complex 
normal form 


Normal forms 
of real subtypes 


Equivalences 


Values of k 






{A+) 
-x'^+i {AD 


At 

for even k 


k > 1 


Dk 








k>4 


Ee 


x^ + y^ 


x^+y^ iE+) 
x^-y^ [E^ 






Er 


x^ + xy^ 


x^ + xy^ 






Es 


x'^ + y^ 


x^ + y^ 







Using the Singular library classify. lib ( Kriiger . 2012[ ) for the complex 
classification and the one-to-one correspondence between the main real singularity 
types and the complex types, the algorithmic classification of a real germ boils down 
to determining to which of the corresponding subtypes the germ is equivalent. For 
the singularity types E-j and E^, there is nothing left to do because each of these 
types has only one real subtype. The rest of the cases is considered one by one in 
the following subsections. 

Throughout the rest of this article we write / for the given input polynomial, 
g for its residual part which can be obtained by applying the Splitting Lemma, and 
c for the corank of /. We also assume that /, and thus g, is a polynomial over Q. 
With these notations, g is a polynomial in c variables. 

4.1. Ax. If c = 0, then it follows that / is of type A^ or of type A^ depending 
on the inertia index A of /. If the inertia index of / is nonzero and less than the 
number of variables in the base ring, then / is both of type A^ and A^ , depending 
on how one chooses to order the variables. We treat those cases as Af . If the 
inertia index is equal to the number of variables in the base ring, / is of type A^ . 
Lastly, if the inertia index is 0, then / is of type A^ . 

4.2. Ak, fe > 1. If c = 1, then the singularity is of type At or of type Aj^ for 
some fc > 1. Furthermore g is a univariate polynomial in this case, say g G Q[x]. 
The value of k is given by the order of g minus 1. This follows since zLx'''^^ and g 
are right-equivalent and thus have the same order. 

Note that if k is even, then At ~ A^^ . Again, we consider those cases as At 
for the sake of simplicity. Let k be odd. Then the sign of the singularity type is 
determined by the sign of the coefficient of x'^'^^ . This follows since it follows from 
Proposition El that jet(5,fc-f 1) = ±(</)i(a;))'=+i = ±{ax)''+^ , where (l){±x''+^) = 5, 
a € M, and the sign depends on the singularity type. Since fc -I- 1 is even and a G M, 
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<j) does not change the sign of the coefficient of x'^^^. We use Algorithm [3l after 
applying the Splitting Lemma in case c = or c = 1. 



Algorithm 3 Algorithm for the case Ak 

Input: / G Q[a;i, . . . , a;„] of complex singularity type Ak, the output polynomial g 

after applying Algorithm [21 the corank c of f and the inertia index A of / 
Output: the real singularity type of /. i.e. A'^ or yl^T, fc € N 



1: if c = then 

2: if A < n then 

3: type := A+ 

4: else 

5: type := A^ 
6: if c = 1 then 

7: k :— ord(5) — 1 

8: if k is even then 

9: type := A+ 

10: else 

11: s :— coefficient of x'^ 

12: if s > then 

13: type A+ 

14: else 

15: type := A~ 

16: return type 



For the rest of the paper we turn our attention to singularities of corank 2. In 
these cases 7^ 5 S m'^ is a polynomial in two variables, say g e Q[x,t/]. Using 
the Singular library classify. lib, we determine the complex singularity type 
and thus the main real singularity type of g, or equivalently /. The purpose of the 
remaining algorithms in the paper is to determine the correct real subtype of g, 
or equivalently /. We now consider each complex type, or equivalently every real 
main type, seperately. 

4.3. D4. The normal form of the complex singularity type D4 is x^y + y^, 
which splits into x^y + y^ {D^) and x^y — (D^) in the real case. The two cases 
can be distinguished by factorization; the details are carried out in Algorithm U) 
Since the determinacy of D4 is 3, it suffices to look at the 3-jet. The number of 
factors of the 3-jet is an invariant of the real subtype which is 1 in the case and 
3 for D^. 

However, using the Singular command f actorize in order to determine the 
number of factors is problematic because the factorization over R differs from those 
over Q and C in some cases. As an alternative, we dehomogenize the 3-jet and count 
the number of real roots of the resulting univariate polynomial which is exactly the 
same as the number of factors of the 3-jet over M. 

If we want to dehomogenize the 3-jet via x t-^ x, y t-^ 1 without reducing its 
degree, we first have to make sure that the coefficient of x^ is non-zero. It is easy to 
check that this is achieved by lines 2 to 13 of Alg orithm [H Fo r the implementation 



in Singular, we used the library rootsur.lib ( Tobid . 2012 ) to count the number 
of real roots of a univariate polynomial. 
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Algorithm 4 Algorithm for the case D4 

Input: (7 g m'^ C Q[a;,y] of complex singularity type D4 
Output: the real singularity type of g, i.e. or 

1: h := jet(.g,3) 

2: Si :— coefficient of in h 

3: S2 coefficient of in h 

4: if (si = 0) then 

5: if (s2 ^ 0) then 

6: swap the variables x and y in h 

7: else 

8: ti coefficient of x^y in h 

9: t2 := coefficient of xy^ in h 

10: if ih +t2T^0) then 

11: apply X 1-^ X, y 1-^ X + y to h 

12: else 

13: apply X X, y 1-^ 2x + y to h 

14: apply X i-^ X, y i-^ 1 to h 

15: n := number of real roots of h 

16: if (n < 3) then 

17: return 

18: else 

19: return 



Remark f 0. Geometrically, the dehomogenization in Algorithm |4] corresponds 
to blowing the 3-jet up at the origin plus choosing a chart. Since the 3-jet is 
homogeneous, blowing-up always yields three lines in the complex case. In the real 
case, however, we get either one or three lines depending on their position w.r.t. the 
real subspace in the complex picture. All the lines lie in the chosen chart because 
the coefBcient of x^ is non-zero. 

4.4. Dk,k > 4. For the cases Dk with fc > 4, the complex normal form is 
x^y + y^~^- It splits up into x'^y + y''~^ (D^) and x'^y — y^~^ (-P ^) for each k 
over the reals. We use the following two results from lSiersmal ( 19741 ) to distinguish 
between the two cases: 

Lemma 11. A singularity of type or is [k — l)-determined. 

Lemma 12. Let j > 4. Then 

x^y -t- a^x^ + aix^^^y + • • • -|- Ojy^ ^ x^y + ajy\ oq, . . . ,aj G M, 

using the R-algebra automorphisms 

X X + pi, where pi = —-{aix^~'^ + ■■■ + aj-iy-'~'^) , 

y y +P2, where p2 = -aox^~'^ . 

By Lemma[TT]the determinacy of a singularity of main type Dk is fc — 1. There- 
fore we only need to consider the (fc — l)-jet of g in this case. Using Proposition [S] 
and Lemma [7j we can transform g into a polynomial of the form 

x'^y + terms of degree higher than 3 
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by factorizing the 3-jet of g as 3^32, 51 and 32 of degree 1, and then applying the 
automorphism defined by f;i 1— )> a;, 52 >— > 2/ to g. We now systematically consider the 
terms of each degree 3 < j < k. By applying the transformations in Lemma [T^ for 
each j, the only term of total degree j which possibly remains is ajyK This term 
vanishes for j < k — \ and it does not vanish for j = fc — 1, otherwise g is not of 
complex type D^. Thus, after applying these transformations, we can write g as 
g = x^y + ay'^~^ with a ^ 0. Clearly if a > then x^y + ay'^~^ ~ x^y + y'^~^ and 
if a < then x^y + ay'^^^ ^ x'^y — y*^"^. 



Algorithm 5 Algorithm for the case Dk, fc > 4 

Input: (7 e m'^ C Q[a:,j/] of complex singularity type Dk, fc G N, fc > 4 
Output: the real singularity type of g, i.e. or 

1: fc := n{g) 

2: h := jet(5,fc - 1) 

3: factorize jet(/i, 3) as h\h2, where hi and /i2 are linear 

4: apply hi ^ X, h2 ^ y to h 

5: for (j = 4, . . . , fc - 1) do 

6: if (jet(/i, - x^y ^ 0) then 

7: write jet(/i, j) — x^y as aoa:-' + aix^~^y + • • • + ajy^ , clq, . . . aj 6 Q 

8: apply a; H> a; — i(aix^^^ + • ■ • + aj^iy^^'^), y ^ y — apa;^^^ to ft, 

9: /i := jet(/i, fc - 1) 

10: write h as h — x^y + ay^~^, / a G Q 

11: if (a > 0) then 

12: return Dj!^ 

13: else 

14: return 



4.5. £^6- In this case, whose complex normal form is x^ + y'*, we have that 
either g^x^ + y'^ {^t) or g ~ a:"^ — y"' {Eq). Therefore there exists an K-algebra 
automorphism (f) of E[a;, y] such that 0(g) = ((/)(a:))^ + (^(y))'* or such that 4>{g) = 
(0(a;))'^ — {(t>{y)Y . Since the coefficients of x^ and y^ in g cannot both be zero, 
we can ensure that the coefficient of x^ is non-zero by swapping the variables if 
necessary. Now, using Proposition [6] and Lemma [71 jet(g,3) factorizes as c{gif', 
c G Q and gi — b^x + biy G Q[x, y]. Again, using Proposition [6l it follows that by 
applying x i— >■ ^~ , y 1-^ y to g, (j) is transformed such that (f>i{x) — c'x, c' G M. 
Since (j) is an automorphism, 0i(y) = cIqx + diy, do,di G M, with di ^ 0. Hence 

{(j){y)y — d\y'^ + terms of degree 4 and higher, not of the form ay*, a G M. 

If we can show that {(f){x)Y does not contain a term of the form ay**, a G M, then 
we can determine whether g is of type or E'^ by considering the sign of the 
coefficient of the monomial y*. Now 

jet((0(a:))3,4)-jet((0(x))3,3) = 3(0i(a;)2)(02(a;) - 0i(a;)) 

= 3(c'a;)2((^2(a:)-0i(a:)), 



which means that {4){x))'^ does not have terms of the form ay"', a G M 
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Algorithm 6 Algorithm for the case Eq 

Input: g Gm^ C Q[x, y] of complex singularity type Eq 
Output: the real singularity type of g, i.e. -Eg" or Eq 

1: h := jet(3,3) 

2: s := coefficient of in h 

3: if (s = 0) then 

4: swap the variables x and y 

5: factorize h into linear factors over Q[x, y], with a factor gi = box + b\y 
6: apply X ^ y^ytog 
7: d := coefficient of y^ in g 
8: if {d > 0) then 
9: return E^ 
10: else 

11: return Eq 
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